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ABSTRACT: A very weakly coupled scalar field with mass m and initial vac- 
uum expectation value V will provide enough mass to close the universe provided 
V — (3 x 10^ GeV)(100 GeV/m) 1 / 4 . We discuss possible models in which such a 
field could arise. 
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Particle physics has provided a number of plausible dark matter candidates, most no- 
tably massive but light neutrinos [1], the lightest supersymmetric particle [2-4], and the 
axion [5-7]. These last two would be cold dark matter, the most plausible type. (For a 
review of dark matter, see Ref. [8].) The axion is especially interesting; a closure density 
of axions does not arise from thermal processes (as would be expected in a hot early uni- 
verse), but rather from a coherent state of particles which approximate a classical field 
with a definite value. Thermal effects do play an important role, however. At very high 
temperatures, the axion field has no potential energy, and so assumes an arbitrary initial 
vacuum expectation value (VEV) V, determined by the details of the initial conditions. 
In most models of early cosmology, this arbitrary value is essentially random and unpre- 
dictable; all we know is that it lies on a circle of circumference 2nf, where / is the axion 
decay constant, related to the axion mass m a via m a ~ (100 MeV) 2 //. As the temperature 
drops, however, a potential turns on (proportional to T~ 8 ), and eventually the axion field 
begins to oscillate about the minimum determined by the potential. The contribution Q of 
these coherent axions to the mass density of the universe (in units of the critical density) 
is given by [5-7] 

nh^o.sf * V( m « ), (i) 

V10 12 GeVy VlO- 5 eV/ V ; 

where h is the current value of the Hubble parameter in units of 100 km/s-Mpc. (For a 
more detailed analysis, see Ref. [9].) Because the initial axion VEV V lies on a circle of 
radius /, it is likely that V ~ /. Putting this into eq.(l) and using the relation between / 
and m a gives an upper limit on / of about 10 12 GeV. 

In this paper we examine the case of a scalar field with mass m and VEV V which we 
imagine to be so weakly coupled that thermal effects can be neglected entirely. We find 
that the contribution of such a field to the mass density of the universe is 

m ^ ' 7 (^IPG^) 2 (lO^W) 1/2 ' (2) 

Note that O ~ V m ' , in contrast to the V m dependence in the case of the axion. While 
this formula is easy to derive, to our knowledge it has never appeared in the literature. 

Are such fields likely to occur? We note that string theory (with four uncompacti- 
fied dimensions) always has at least four massless scalar particles (the string dilaton, the 
breathing mode of the compact dimensions, and two partners demanded by supersymme- 
try) [10]. These particles can get mass only after supersymmetry is broken. It is widely 
anticipated that the spontaneous breakdown of SU(2)xU(l) will turn out to be related 
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to supersymmetry breaking, yielding natural masses of order the weak interaction scale, 
~ 100 GeV, for these scalars. However, there is no reason to believe that these scalars 
will emerge from the early universe (after a period of inflation, say) at the minima of the 
potentials generated by supersymmetry breaking. They may have large, arbitrary VEVs. 

Indeed, a very similar proposal has been made by Linde [11], but to provide inflation 
rather than dark matter. Linde envisions a region of space in which a scalar field with 
mass m happens to have a large VEV V, large enough that the stored energy density 
dominates over that of all other particles and fields. This energy density drives inflation. 
Linde refers to this as "chaotic inflation," since the field takes on a chaotic pattern of 
values throughout the universe. We envision a similar scenario, except that we assume 
the energy density stored in our field happens to be much less than the thermal energy 
of ambient hot particles, which we imagine were produced by reheating after a period of 
inflation. Our field does not drive inflation, but it does provide for dark matter. Given 
the similarity to chaotic inflation, we refer to this scenario as chaotic dark matter. 

Of course, to be a dark matter candidate our particles must be stable or nearly so, 
with a decay rate no larger than the Hubble parameter. Also, they must be self-coupled 
weakly enough to avoid efficient transfer of energy out of the coherent state. As we will 
see, this latter constraint is not a problem, but the former is an impediment to identifying 
our field as the dilaton or one of its cousins. We would generically expect a dilaton with 
mass m to have a decay rate of order m 3 /mp l7 where m P i = 1.2 x 10 19 GeV is the Planck 
mass. (This is because the dilaton should couple with gravitational strength to all other 
fields.) Requiring this to be less than the Hubble parameter translates into m < 100 MeV, 
a factor of a thousand smaller than the weak interaction scale. We will not attempt to 
solve this dilemma here; like the field which is hypothesized to drive chaotic inflation, our 
field need not be one of the string fields mentioned above. Producing a closure density 
with a mass of 100 MeV requires a VEV of 2 x 10 9 GeV. 

We now turn to demonstrating eq.(2). For the most part we can follow the analysis 
originally used for the axion [5-7]. The equation of motion for our field <p is 

(p - R~ 2 V 2 v + 3(R/R)<p + m 2 ip = 0, (3) 

where R(t) is the scale factor; R ~ t 1 / 2 when the universe is radiation dominated and 
R ~ when the universe is matter dominated. We take the initial value of <p(x, t) to 
be V everywhere; possible spatial inhomogeneities will be treated later. Until the age t of 
the universe is comparable to the period 1/m of the oscillations, the field <p will not have 
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had enough time to change significantly, and hence remains frozen at <p(t) = V. Once 
t becomes of order 1/m, however, (p(t) begins to oscillate. Call the temperature at this 
time T\. Afterward, the time-averaged energy density in the field is p = \m 2 ip 2 , with 
(f{t) = V when T = T\. This energy density decreases like the energy density of particles 
with zero momentum; that is, it decreases as l/i?(t) 3 , or equivalently as T 3 , where T is 
the temperature. We thus find 

p = lmV 2 (^m) 3 , (4) 

where Tq = 2.75 K is the present temperature. We compute T\ from the usual time- 
temperature relation, t = m P i/(7T 2 ), where 7 = 34.3 in the standard model and 7 = 49.4 
in the supersymmetric standard model. Using t = 1/m and 7 = 34.3 we find 

T -< 6xl ° 9GeV >(io(nW) 1/2 ' < 5 » 

Putting this into eq.(4) and dividing by the critical density p c = 3H 2 /8nG yields eq.(2), 
our central result. 

We shall now consider self interactions of the (p field. In keeping with the idea that its 
potential is generated by supersymmetry breaking, we parameterize the interactions as 

00 

Ant(^) = -m 2 Xn ™p/ n 2 V n (6) 

where the X n are dimensionless numbers, assumed to be of order one. (Derivative inter- 
actions could also be included, and would lead to similar results.) These interactions will 
allow energy to be transferred out of the zero-momentum particles and into particles with 
relativistic energies; the energy will then red-shift away faster, like T 4 instead of T 3 , for 
as long as the particles remain relativistic. 

Following the analysis of Ref. [6] for the axion, we consider a perturbation of the form 

(p(x, t) = A{t) cos(mt) + A k (t)e tls (7) 

where Af, is treated as infinitesimal. Substituting in eq.(3) and taking R and A to be 
approximately constant, we get 

Ak(t) + [(k/R) 2 + m 2 + 6X 3 m 2 (A/m Pl ) cos(mt) + . . ]A k (t) = 0, (8) 

where the ellipses stand for the higher order terms in eq.(6). These will be suppressed 
by additional powers of A/m Ph and can be neglected. Eq.(8) is Mathieu's equation. For 

4 



certain bands of values of k, Ay. grows exponentially. This would mean that there is an 
instability which could pump energy into high momentum modes. However, as noted 
above, we need only be concerned with modes for which k m, and these modes account 
for an exponentially small range of values of k at large k. Thus the unstable modes do not 
have enough time to grow significantly before they are redshifted into the regime of stable 
wavenumbers. Therefore, just as for the axion, this effect can be ignored. 

To conclude, we have shown that a coherent scalar field with mass m and initial 
VEV V contributes to the mass density of the universe according to eq.(2). A field with 
m ~ 100 GeV and V ~ 3 x 10^ GeV would provide a critical density of cold dark matter. 
Such a field might arise is certain string models, but generically the corresponding particles 
would have lifetimes that are too short unless the mass is reduced to m < 100 MeV. 

This work was supported in part by NSF Grant No. PHY-86-14185. 
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